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ABSTRACT 

We examine the idea that diffuse and giant molecular clouds and their substructure 
form as density fluctuations induced by large scale interstellar turbulence. We do this by 
Q\ I closely investigating the topology of the velocity, density and magnetic fields within and 

■ at the boundaries of the clouds emerging in high-resolution two-dimensional simulations 

^ ' of the ISM including self-gravity, magnetic fields, parameterized heating and cooling 
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O ! and a simple model for star formation. We find that the velocity field is continuous 

across cloud boundaries for a hierarchy of clouds of progressively smaller sizes. Cloud 
boundaries defined by a density-threshold criterion are found to be quite arbitrary, with 
P_|! no correspondence to any actual physical boundary, such as a density discontinuity. 

Abrupt velocity jumps are coincident with the density maxima, indicating that the 
^ ■ clouds are formed by colliding gas streams. This conclusion is also supported by the 

c/3 . fact that the volume and surface kinetic terms in the Eulerian Virial Theorem for a 

cloud ensemble are comparable in general. The topology of the magnetic field is also 
suggestive of the same process, exhibiting bends and reversals where the gas streams 
^ ' collide. However, no unique trend of alignment between density and magnetic features 

■ is observed. Both sub- and super- Alfvenic motions are observed within the clouds in the 

simulations. 

In the light of these results, we argue that thermal pressure equilibrium is irrelevant 
for cloud confinement in a turbulent medium, since inertial motions can still distort or 
disrupt a cloud, unless it is strongly gravitationally bound. Turbulent pressure con- 
finement appears self-defeating, because turbulence contains large-scale motions which 
necessarily distort Lagrangian cloud boundaries, or equivalently cause flux through Eu- 
lerian boundaries. 

We then discuss the compatibility of the present scenario with observational data. 
We find that density-weighted velocity histograms are consistent with observational line 
profiles of comparable spatial and velocity resolution, exhibiting similar FWHMs and 
similar multi-component structure. An analysis of the regions contributing to each veloc- 
ity interval indicates that the histogram "features" do not come from isolated "clumps" , 
but rather from extended regions throughout a cloud, which often have very different 
total velocity vectors. 
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Finally, we argue that the scenario presented here may be also applicable to small 
scales with larger densities (molecular clouds and their substructure, up to at least 
n ~ lO'^-lO'^ cm— 3), and conjecture that quasi- hydrostatic configurations cannot be 
produced from turbulent fluctuations unless the thermodynamic behavior of the flow 
becomes nearly adiabatic. We demonstrate, using appropriate cooling rates, that this 
will not occur except for very small compressions 10~^ pc) or until protostellar 
densities are reached for collapse. 
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1. Introduction 



Theoretical models of interstellar clouds and clumps 
most frequently assume static or stationary configu- 
rations, the c louds being eit her confined by the exter- 
nal pressure ( [Vlaloney f 988| ; Bertoldi & McKee 1992, 
herafter BM92), or in "virial equilibrium" between 
its self-gravity and some form of internal energy, be 
it thermal (e.g., Chieze 1987), micro-turbulent (e.g., 



Chandrasckhar 1951; 


Bonazzola et al. 1987 




Leorat. 


Passot & Pouquct 1990 




Vazquez-Semadeni 


& 


: Gazol 


1995) or magnetic (e.^ 


r 


Shu, Adams & Lizano 1987; 



ing out of instabilities or coagulation of smaller clouds 

Nevertheless, the 



(see Elmegreen 1993a for a review) 
interstellar medium (ISM) is well known to be highly 
violent (e.g., McCray fc Snow 1979| ), and recently it 



has become increasingly accepted that it is turbulent 



throughout (e.g., Bcalo 1987; see also the volume "In- 



terstellar Turbulence" [Franco & Carramiiiana 1998|) 



In such a medium, clouds mav naturallv form as tur- 



bul bnt density fluctuations as well. 



Within this dynamic, turbulent framework, it is 
important to reconsider the meaning and feasibility 
of pressure cloud confinement, either thermal or tur- 
bulent. This is most adequately done in the context 
of the Virial Theorem (VT), which is a very useful 
tool for describing the balance between all the physi- 
cal ,a£!EiitsJiLJXi£d£iiLLLajLxk^ rnstnTnarily the VT 



is ifnnsidprprl in a T.agra.ngian frame i p a frame 



moving with the flow. However, there are circum- 
stances when one would prefer to use an Eulerian 
description. For example, a Lagrangian description 
is clearly not well suited for describing the propaga- 
tion of a wave, in which the fluid does not move along 
with the perturbation, while in the case of an isolated 
object a Lagrangian description fits most naturally 
(see also Vazquez- Semadeni, Passot & Pouquet 1996, 
hereafter VSPP96). In the interstellar case, if signif- 
icant mass exchange occurs between a cloud and its 
surroundings, such as an accreting cloud, an Eulerian 
(fixed in space) may be preferable since a Lagrangian 
boundary would undergo severe distortions. Parker 
(1969) derived an Eulerian form of the Virial Theo- 
rem (EVT) in tensor form, but he explicitly neglected 
the mass flux through the surface of the cloud. More 
recently, McKee & Zweibel (1992) (hereafter [MZ92| ) 
have written all the terms entering the EVT: 



-Ie — 2 I £th + £ kin 



%.h — Tv\ 



kin 
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1 d$ 
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where Ie = jy pr^dV is the moment of inertia of the 
cloud, fth = 3/2 Pt\idV is the thermal energy, with 
Pth being the thermal pressure; ^kin = 1/2 Jy pu^dV 
is the kinetic energy, T^h = XiPthfiidS is the 

surface thermal term, Tkin = 1/2 §gXipuiUjhjdS 
is the surface kinetic term, M = I/Stt JyB^dV + 
Jg XiTijfijdS is the magnetic term, with Tij the Maxwell 
stress tensor, W = Jy pXiQidV is the gravitational 
term and <^ = §g px'^UihidS is the flux of moment of 
inertia through the surface of the cloud. Sums over 
repeated indices are assumed unless otherwise stated. 
Frequently, the surface terms are neglected alto- 



gether, especially in observational work (Larson 1981 



Torrelleset al, 1983; Myers fc Goodman, 1988a; Fuller 



Myers 1992), although it is also a common practice 



in theoretical studies, as a consequence of assuming 



isolated clouds (e.g., Chandrasckhar & Fermi 1953 



Parker 1969, Parker 1979). The most notable excep- 
tion is the thermal pressure surface term 7^h, which 
is frequently invoked for "pressure confinement" (e.g. 



McCrea 1957; Keto fc Myers, 1986; Maloney 1988 



BM92; McLaughlin fc Pudritz 1996; Yonekura et al 



1997). By analogy, MZ92 have considered the possi- 



bility of turbulent pressure confinement by means of 
the term T^in- However, in the present paper we will 
argue that both kinds of pressure confinement require 
idealized conditions that are not likely to be realized 
in the actual ISM. If the interstellar medium (ISM) 
is globally turbulent, the velocity field is in general 
locally nonzero, and thermal pressure balance may 
be irrelevant, since a cloud can still undergo defor- 
mation due to the inertial motions, which will induce 
a nonzero RHS via the kinetic terms i?kin, "^kin, and 
d^/dt. 

Concerning turbulent pressure confinement, its fea- 
sibility requires the underlying assumption that the 
turbulence be microscopic, so that it can be consid- 
ered isotropic compared to the scale of the cloud. 
However, turbulence is inherently a multi-scale phe- 
nomenon, and is expected to contain excitation at 
scales comparable or even larger to that of the cloud. 
It is in this context that the clouds may be considered 
as density fluctuations produced by larger-scale com- 



pressive turbulent motions (e.g.. Hunter 1979; Lar- 
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son 1981; [Hunter fc Fleck 1982| ; [Hunter et al. 198e| ; 
Elmegreen 1993b; Vazquez-Semadeni, Passot & Pou- 



quet 1995, hereafter VSPP95). Such clouds are then 



turbulent fluctuations, and can either rebound, frag- 
ment or collapse depending on the compressive energy 
available, the cooling ability of the flow, the topology 
of the velocity and magnetic fields, the production of 
internal turbulent motions and a number of instabil- 
ity mechanisms (e.g., Elmegreen fc Elmegreen 197^ ; 
Vishniac 1983, 1994; [Hunter et al. 1986[; Tohline, Bo- 
denheimer fc Chr i stodoulou 1987; Stevens ct al. 199^ ; 



Elmegreen 1993b; VSPP96; Kornrcich & Scalo 1 



Note that in this scenario, clouds and clumps need 
not be in a static or quasi-static equilibrium at any 
time. Moreover, the fact that clouds and clumps are 
observed to be internally turbulent suggests that the 
turbulence has not been completely dissipated inside. 
In fact, it is known that shocks running obliquely 
through a densi ty gradient will produce further inter- 
nal turbulence ( Klein et al. 1995| ; Kornrcich fc Scalo 
19c(|). 

In this paper, we analyze the density, velocity and 
magnetic fields in two-dimensional numerical simu- 
lations of turbulence in the ISM in order to estimate 
the importance of the transport processes through the 
cloud boundaries. We start in §|| by discussing the 
role of the kinetic terms in the Virial Theorem and 
the implications for pressure confinement. Then, after 
briefly describing the numerical simulations (§^), we 
proceed to a description of our results (§ ^), starting 
with a discussion of the density, velocity and mag- 
netic features that arise, and their spatial correla- 
tions (§4.1). Then we present an evaluation of the 
surface and volume kinetic terms in the virial theo- 
rem, suggesting that their similarity may be indica- 
tive that both types of terms are representative of the 



same phenomenon (§ 4.2). The super- or sub-Alfvenic 
character of the motions is discussed in §4.3. In § [5[ 



we present several comparisons with known observa- 
tional data, such as line profiles observed in cloud 
co mpl exes (§5J^), the topology of the mag net ic field 
(§ and the lifetimes of the clouds (§ x3). In § 

H we then discuss several issues and possible caveats, 
like whether the scenario proposed here can be ap- 
plicable to smaller, denser scales, the low probability 
of forming quasistatic clumps in a turbulent medium, 
and the effect of the dimensionality of our simulations. 
Finally, we give a summary and some conclusions in 



be regarded as theoretical justification for the sug- 
gestion made over 45 years ago by Chandrasckhar & 



Munch (1952) that an alternative to the picture that 



visualizes the interstellar medium as consisting of a 
distribution of discrete clouds might be necessary. 

2. Kinetic Terms in the Eulerian Virial The- 
orem 



As is well known (e.g., [Parker 1969| ; [MZ92[ ), the 
EVT can be obtained from the conservative form of 
the momentum equation. 



d{pui) dipu^Uj) 



dt 



dxi 



(2) 



(where F}^ generically represents the various force den- 
sities involved in the problem) by dotting it with 
the position vector x and integrating over a volume 
V fixed in space. We remind the reader that the 
full Virial Theorem, eq. (|l|), holds always because 
it is a direct consequence of the momentum equa- 
tion, although this does not necessarily imply that 
the system is in equilibrium or, in particular, that the 
linewidth of an interstellar cloud is directly related to 
its mass and radius ( Maloney 1988[ ). A detailed anal- 
ysis of the VT in our simulations will be presented 



elsewhere (Ballesteros-Paredes et al. 199S 



In the present paper we focus on the relative im- 
portance of the surface and volumetric kinetic terms, 
and their role in shaping the clouds. It can be eas- 
ily checked that, in the EVT (eq. Q), the kinetic 
energy 5kin, the kinetic surface term Tkin, and the 
time-derivative terms Ie and d^/dt, originate from 
the LHS in the momentum equation (eq. 1^). While 
the meaning of the kinetic energy term £'kin is ob- 
vious, the kinetic surface term Tkin has two possible 
interpretations. One is as half the flux of the instan- 
taneous rate of change of the moment of inertia pXiUi 
through the cloud's surface. Alternatively, it can be 
interpreted as the sum of the ram pressure plus the 
kinetic stresses, both evaluated at the surface of the 
cloud: 



XipUiUjUjdS — (p Xipu fiidS -f 
s J s 



XipuiUjhjdS, 



(3) 



It is interesting to note that the present paper may 



where i ^ j in the last surface integral. The first term 
on the RHS of the eq. (||) is analogous to the thermal 
pressure surface term. The second term, on the other 
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hand, reflects the fact that the turbulent motions are 
not isotropic, giving off-diagonal contributions to the 
total pressure tensor Ily = Pth^ij + pUiUj (Landau 
& jifshitz 1987). This exhibits in a clear way the 



difference between the isotropic nature of the thermal 
pressure and the anisotropic nature of the "turbulent" 
pressure. 

Another important difference between the thermal 
and "turbulent" pressures lies in the fact that the 
latter is in general expected to involve motions of 
scales comparable to that of the cloud, as mentioned 
in But since the velocity field carries mass with it, 
such large-scale motions necessarily imply a mass flux 
across fixed, Eulerian cloud boundaries. Conversely, 
in a Lagrangian description, such turbulent motions 
should lead to severe distortions of the Lagrangian 
cloud boundary, possibly giving it a fractal characte r 
(e.g., Scalo 1990; [Falgarone, Phillips fc Walker 199l|) . 

The exception to the scenario depicted above is the 
case of strongly self-gravitating clumps in (magneto- 
) hydrostatic equilibrium, from which the gas can- 
not escape once it has been captured ("decouples" 
from the intercloud medium). In this case, the cloud 
should be bounded by an accretion shock and, al- 
though strictly speaking there is flux across the shock, 
the accreted material is not mixed far beyond the 
shock into the cloud's "body". However, it seems 
to us that the question as to whether such quasi- 
hydrostatic configurations can actually be produced 
in a turbulent ISM remains open, as discussed in § 

o. 



3. Numerical Simulations 

In the following sections, we discuss the topology 
of the density, velocity and magnetic fields in 2D nu- 
merical simulations of the ISM based on the model of 
Pass ot, Vazqu ez- Semadeni & Pouquet (1995) (here- 
after |PVSP95D , which represents 1 kpc^ of the ISM 
on the Galactic plane, centered at the solar Galac- 
tocentric distance. We refer the reader to that pa- 
per for the equations and parameters of the model, 
which includes self-gravity, magnetic fields, parame- 
terized cooling and diffuse heating, the Coriolis force, 
large-scale shear, and parameterized localized stellar 
energy input due to ionization heating. The param- 
eterized cooling is as in Chiang & Bregman (1988)1 , 
who fitted piecewise power laws to the standard cool- 
ing calculations of Dalgarno fc McCray ("l972] and 



VSPP96, the net result of the combined cooling and 
diffuse heating laws is that the flow effectively be- 
haves as a polytrope, with an effective polytropic ex- 
ponent 7cff which is piecewise constant over the tem- 
perature range spanned by the flow. Alternatively, 
it was shown in VSPP95 that, for the heating and 



cooling functions used there, the thermal equilibrium 
values of the temperature and pressure can be ex- 
pressed as functions of the density in the range 100 K 
< T < 10^ K for non star- forming regions. Therefore, 
7off can also be expressed as a piecewise constant over 
the corresponding density r ange. For the fiducial val- 
ues ad opted in the model ( PVSP95 ), 7cff takes the 
values (|VSPP96|) 



7off 




1.57 < p (100 < T < 2000) 
0.39 < p< 1.57 (2000 < T < 8000) 
3.15 X 10-3 < p < 0.39 (8000 < T < 10^) 
4.28 X 10-2 < p < 3.15 X 10-3 
(10^ < T < 4 X 10^), 

(4) 

where densities are in units of cm ^ and temperatures 
in Kelvins. In eq. (|^), equilibrium temperature ranges 
equivalent to the specified density ranges are indi- 
cated in parentheses. In the simulations, stellar ion- 
ization heating is modeled by means of a local heat- 
ing source which is turned on if the density exceeds a 
critical value pc = 30 (p), where the brackets denote 
an average over the whole volume of the simulation. 
Note that temperatures over a few x 10* K are never 
reached in the simulations, because supernovae are 
not included, since these require non-trivial modifica- 
tions to the code ( Gazol fc Passot 1998 ). The max- 
imum densities reached in the simulations are ~ 100 
cm-3, with T - 100 K. 

The simulations we use are refinements of the 



PVSP95| model at a higher resolution (800x800 grid 



points), presented in Vazquez-Semadeni, Ballesteros- 
Paredes, & Rodriguez (1997), hereafter |VSBPR97 . 
The stellar energy injection maintains the turbulence 
in the medium, promoting further cloud formation. 
Nevertheless, we turn off star formation shortly be- 
fore the time at which the data are analyzed, in or- 
der to allow for the largest possible density gradients, 
since otherwise the stellar heating prevents the den- 
sity from reaching values significantly larger than pc 
by causing the neighboring gas to re-expand. The var- 
ious physical quantities are in units of po = 1 cm-^, 
Wo = 11-7 km s-\ To = 10* K and = bpG (see 



Raymond, Cox & Smith (1976). As discussed in 



PVSP95) 



We have developed a cloud-finding numerical algo- 
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rithm in physical space (as opposed to the position- 
velocity space of observational data) for the density 
data from the simulations ( pallcstcros-Paredes and 
Va: :quez-Semadeni 1995| ). A cloud is defined as a 



connected set of points (pixels) whose densities are 
larger than a given threshold px- The clouds obtained 
through this procedure are clearly quite arbitrary, 
since as is increased the boundary of the cloud 
simply moves "inwards" . However, eventually an in- 
crease in pt will lead to a more qualitative change 
in which a given cloud splits into two or more "chil- 
dren", in a similar manner to the "structure trees" 



used by [Houlahan fc Scalo (1992)| . This procedure is 
in a sense analogous to performing observations using 
different tracers sensitive to different density ranges. 
Alternative cloud identification algorithms based on 
laocating density maxima (e.g., Williams, de Geus & 
Blr :z 1994 ) were not investigated, but we do not ex- 
pect that such procedures would affect the qualitative 
nature of our results, which, as will be seen, appear to 
follow quite generally from the physics of the problem. 

Given the above algorithm for defining the clouds, 
the kinetic terms entering the VT are then evaluated 
numerically for each cloud, having previously sub- 
stractcd its bulk mass-averaged velocity, defined as 



(Ui)r 



IvPdV ' 



(5) 



in order to measure only the contribution associated 
to the fluctuations. We refer to this as measuring the 
velocities in the "cloud frame" . 

4. Results and Discussion 

4.1. Density and velocity field topology 

In figs. |l|a, b and c we show three hierarchical levels 
of the density field (gray scale) , with the correspond- 
ing velocity fields (arrows). In fig. 0c we further de- 
note super- and sub-Alfvenic velocities with black and 
white arrows, respectively (see § where the Alfven 
speed is defined as v"^ = B"^ / p in the code's units. 

Figure |l|a shows a snapshot of the simulation at 
t — 7.2x 10''yr. The density field ranges from 0.01 
to ~ lOOcm^'^. Figs. |^b and c show two subsequent 
hierarchical levels of clouds, obtained a.i p^ — A and 
8, respectively. The sizes of the boxes shown in the 
three panels are 800, 215 and 45 pixels per side (1 
pixel = 1.25 pc). We refer to the cloud in fig. |l|b as 
the "parent" cloud, and to that in fig.[|c as the "child" 



cloud. Note that only fig. |i|c shows the velocity field 
in the cloud frame. 

Inspecting the plotted density and velocity fields, 
we observe that the fluid velocity at the cloud bound- 
aries is in general continuous, at any level of the hi- 
erarchy. Sharp changes in the velocity field (shocks)]^ 
are in general oblique and tend to occur mostly at 
the centermost parts of the clouds, where the highly 
filamentary density field exhibits "ridges" . This indi- 
cates that the density features are produced by collid- 
ing streams, a fact that can be seen more accurately 
in figs. b and c. These show cuts, along the x- 
(at y — 514, fig. |^) and y- (at x = 602 and x — 614, 
figs, b and c respectively) axes, of the density p (solid 
line), and the x- and y-components of the velocity, Ux 
(dotted line) and Uy (long dashed line) for a region 
containing the clouds of figs. |l]b and c. The axes in 
these plots indicate pixel number. For added clarity, 
we show in fig. the lines along which these cuts are 
taken. 

The central parts (3-5 pixels) of the density peaks 
are seen to correspond in general to large negative 
longitudinal gradients dui/dxi (no sum over i) of 
the velocity, i.e., shocks. However, the clouds (i.e., 
large density values) extend beyond the shocks for 
many more pixels, where the velocity component ex- 
hibits plateaus with superimposed small-scale struc- 
ture. These plateaus are indicative of the approaching 
streams (extended inward motions). As an example, 
consider the peaks labeled "1" and "2" infig.|2^. Peak 
1 corresponds to a shock at a; = 602, and is flanked 
by an expansion wave {dux/dx > 0) on the left, and 
by a plateau with Ua, < on the right, although the 
plateau itself contains substructure. This peak also 
corresponds to a shock along the y-direction (flg. |^b). 
Interestingly, peak "2" in fig. ^ does not seem to cor- 
respond to a shock along the x-direction, but in fig. 
^ it can be seen that it is caused by a shock in the y- 
direction. In order to further illustrate this point, in 
fig. 1^ we show the same cloud as in fig. |^, but show- 
ing the density field as a gray-scale image and the 
divergence dui/dxi of the velocity field with coun- 
tours. Local density maxima clearly coincide with 
local (negative) minima of the divergence. 

It is worth noting that in figs. ||a, b and c, larger- 
scale density features are also seen to correspond to 
large-scale negative gradients of the velocity field, al- 



^In the simulations, shocks are sharp gradien 
3 to 5 pixels, due to the action of viscosity 



jng over 
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though with much substructure. See, for example, the 
large cloud complex (multiple density peak) extend- 
ing from 580 <j x <^ 680 in fig. The very largest- 
scale inward flows are probably of gravitational rather 



than turbulent origin (e.g., Vazquez-Semadeni c 
sot 



Pas- 



1998), because the simulations do not include su- 



pernovae which could induce turbulent motions at 
those scales, but in the real ISM both cases are likely 
to occur. 

There also appears to be a certain amount of spa- 
tial correlation between density and vorticity features. 
In figure ^ we show the vorticity field (contours) su- 
perposed onto the density field (gray-scale) for the 
cloud of fig. I^c. In this case, the correlation between 
local maxima of the density field and maxima (or min- 
ima) of the vorticity field is not as tight, the critical 
points of the vorticity being shifted by a few pixels 
with respect to the density peaks. Nevertheless, there 
is still a clear correspondence between vorticity and 
density features. This correlation can be understood 
because, as mentioned above, the stream collisions are 
in general oblique. It is possible then that the detailed 
structure inside the clouds consists of shocks and tan- 
gential discontinuities. Moreover, the shocks are gen- 
erally curved and are encountering density variations 
in the pre-shock gas, so vorticity production is ex- 
pected behind them (Hayes 1957; Passot & Pouquet 
19^ [Kornreich fc Scalo 1998| ) 

The continuity of the velocity field through the 
cloud boundaries has an interesting subjective impli- 
cation. The definition of a cloud by a density thresh- 
old is actually seen to be rather arbitrary, since the 
clouds do not have sharp physical boundaries. Thus, 
as long as the change in the density threshold does 
not imply a "splitting" of a cloud into its offspring, 
it only amounts to considering the same cloud out to 
different distances from the maximum, but without 
there being any clear "edge" to the cloud. Note that 
this need not be in conflict with observations sug- 
gesting that molecular clouds have sharp edges ( Blitz 
since in this case the observed boundaries may 
just reflect a transition from molecular to atomic gas, 
without a discontinuity in t he gas mass density (e.g.. 



Goldsmith 1987, Blitz 1991) 



4.2. Virial surface terms and pressure con- 
finement 

The velocity flelds observed in the simulations indi- 
cate a strong flux of the relevant dynamical quantities 
through the cloud boundaries at all scales. To quan- 



tify this, in fig. H we show a plot of the kinetic energy 
^kin vs. the kinetic surface term Tkin for the parent 
complex shown in fig. |l|a and all its daughter clouds 
resulting from setting px =4, 8, 11 and 16, and with 
areas larger than 300 pixels. The latter requirement 
is imposed in order to avoid including clouds so small 
that they are strongly influenced by viscous and dif- 
fusive numerical effects. The clouds shown in figs. Eb 
and c are represented with a different symbol in fig. q 
for identification. Note that both kinds of virial terms 
are computed in the cloud frame. 

It is seen that in general the surface kinetic virial 
term is of the same order of magnitude as the volume 
term, i.e., the kinetic energy contained in the cloud's 
volume. This means that both terms are of similar 
importance in shaping and supporting the clouds. At 
a more detailed level, this seems to be a refiection 
of the continuity of the velocity field. The fiow is 
entering the clouds, and shocking at their innermost 
regions. Thus, the two terms refer to essentially the 
same process, only measuring it at different places 
(one over the cloud's volume, the other at the cloud's 
boundary but weighted by the distance to the center 
of mass), explaining their similar values. This is anal- 
ogous to the well known property of the surface and 
volume thermal pressure terms, that if the pressure is 
constant the two terms cancel. 

It is important to remark that the similarity be- 
tween these two terms does not imply that the cloud 
is "confined" by turbulent pressure, in the sense that 
this is not a static configuration. First, the two terms 
do not cancel each other exactly, leaving a net contri- 
bution available for shaping the cloud (the contribu- 
tion £kin ~ ^in to cPl/dt^). Sccond, it is important 
to note that the points in fig. |^ lie sometimes above 
and sometimes below the line £kin = ^in, indicating 
that for some clouds the contribution to cPl/dt^ is 
positive, but negative for others. This implies that 
the kinetic terms sometimes provide net "support" 
and sometimes cause net compression, although in 
general the clouds are expected to suffer distortions 
that cannot be classified as either one of the above. 
Thus, the clouds are evolving with time, continuously 
changing their volume and shape, since they are the 
density fiuctuations produced by the turbulence^ 



^See the accompanying video to VSPP95 for a non 
magnetic example 



An mpeg video of a fully MHD 
miilatinn from |PVSP9,4 among others, can he seen at 



littp: / /www. astroscu.unam.mx/turbulence/movies .html 
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The thermal pressure, as discussed in |VSPP95| , 
shows httle spatial variation (except near star forma- 
tion sites where it is much larger) because < 7eff < 1 
in general for the temperature range spanned by the 
simulations. In particular, 7eff = implies an isobaric 
medium. Note that "feS < 1 is in agreement with the 
elementary fact that in the ISM denser regions are 
colder in that temperature range (Myers 1978). In 



other words, the near constancy of the pressure in 
this temperature range is a consequence of thermal 
balance (giving 7eff < 1) in a medium whose density 
is determined by the turbulent motions. However, 
this near "pressure equilibrium" has no effect in con- 
fining the clouds, w hich are i n a state of constant 
change. As stated in |VSPP95| , the pressure is slaved 
to the density field, which in turn is determined by 
the turbulence. 

4.3. Super- and sub-Alfvenic motions and mag- 
netic field topology 

Another point worth discussing is the sub- or 
super-Alfvenic character of the velocity, since it has 
been traditionally argued that the motions in molecu- 



lar clouds are superson ic but sub-Alfvenic (e.g., 3hu 
Adams & Lizano 1987), although more recently it has 
been claimed by Padoan & Nordlund (1998) that this 
may not be so, but rather that motions may generally 
be super-Alfvenic. Since in our simulations clouds 
form out of the general ISM, they represent a good 
test to see which nf these cnnditinns develop, rather 



thain using some pre-determined a,ssumptimi about 



the sub- or super- Alfvcnic character of the problem 
(e.g., [Mouschovias 1987| ; [Lizano fc Shu 1989| ). In 
fig. ^ the black and white arrows respectively de- 
note super- and sub-Alfvenic velocities. It is seen 
that inside the cloud there are both sub- and super- 
Alfvenic velocities. Since the velocity field does not 
show large magnitude fluctuations at the transition 
sites, this can only be understood in terms of a change 
in the Alfven speed Va. To see this, in fig. |l|c contours 
of the magnitude B of the magnetic field are shown 
as well. The contours span the range from 4 to 10 
liG, in intervals of 1.5 fiG. Thicker lines indicate pro- 
gressively larger values of B. The velocity is seen to 
be sub-Alfvenic everywhere outside the cloud. Inside 
the cloud, it is super-Alfvenic in the lower region of 
the cloud. In the upper region, the velocity is sub- 
Alfvenic due to a combination of large i3-values and 
low densities. Looking down towards the center, an- 
other super-Alfvenic region is seen, due mainly to the 



increase in the density, which reaches p = 55 cm~^ 
at the peak. However, at the peak, the velocity be- 
comes sub-Alfvenic again, this time due to a sharp 
decrease in the velocity magnitude itself - sort of a 
"stagnation point" at the density peak. These vari- 
ations are thus due again to the collisions between 
the magnetized gas streams, which push the field and 
cause magnetic shocks and field reversals where the 
fluid shocks occur, and in general produce a rather 
chaotic flow. One such reversal is seen to occur on 
the filament extending out of the cloud towards the 
lower left corner of fig. |l|c. This is seen on the con- 
tours as a decrease in the field's intensity towards the 
central ridge of the filament. Figure which shows 
the magnetic field vectors on top of the density (gray 
scale) and the vorticity (contours), shows this phe- 
nomenon more clearly. 

Note that this trans- Alfvenic character of the mo- 
tions may explain why the density peaks are "fat" 
compared to the stream-collision sites: the informa- 
tion of the presence of the collision is often able to 
travel upstream of the flow. Furthermore, since the 
collisions are generally oblique, it is likely that magne- 
tosonic waves propagating perpendicular to the col- 
lision surface may overtake the perpendicular com- 
ponent of the oblique fluid motion, even if the total 
fluid velocity is super-Alfvenic. In contrast, it is well 
known that in Burgers turbulence, where no pressure 
is available, the density structures are as thin as al- 



lowed by viscosity exclusively (see, e.g., Chappell & 
Scalo 1998| ). 



The advection of the magnetic fleld described above 
is contrary to the usual assumption that the fleld con- 
trols the fluid motions, forcing the gas to flow along 
the fleld lines. Nevertheless, it is noteworthy that 
the field strengths developed by the simulations in 
both the intercloud medium and in the clouds are 
not unrealistic. The lowest values of the field oc- 
cur at the intercloud medium, and are ~ 1/^G, while 
the largest values occur in the densest regions, with 
values ^ 25fiG. Nevertheless, the field is quite inter- 
mittent, and does not follow a unique scaling with 
the density ( PVSP95 ), exhibiting significant varia- 
tions within the cloud of fig. Moreover, even un- 
der near-equipartition conditions between kinetic and 
magnetic energy, the velocity field should be expected 
to exert as much of an influence on the magnetic fleld 
as viceversa. In § 5.2 we discuss the degree to which 
this fleld topology is consistent with observations. 
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5. Comparison with observations. 

The scenario presented here differs significantly 
from the standard view in which clouds and clumps 
have well defined boundaries (density discontinuities) 
which separate them from a much more tenuous, 
warmer me dium, as in simulations o f cloud coUi- 
sions (e.g., Pumphrey fc Scalo 1983| ; Klein ct al 
19E|5|; Miniati et al. 1997 ), of cloud con finement (e.g. 



BM92 ), quasi-static cont raction (e.g., Lizano fc Shu 
19^91) and stabihty (e.g., [Ebbert 19"55| ; |Bonnor 195e| ; 
McLaughhn fc Pudritz 1996| ), etc. On the other 
hand, Scalo (1990) has suggested that such a standard 
model may be an over-idealization stemming from the 
low resolutions and column density ranges of early 
studies, and proposed instead a much more complex 
scenario in which clouds do not have a well-defined 
identity. Our results support this view, and in fact 
such poorly-defined identity is reflected in the arbi- 
trariness of cloud- and clump-finding algorithms (e.g., 
WiUiams, de Geus fc Blitz"T99^ ; |VSBPR97D . It is nec- 
essary, however, to test whether the scenario proposed 
here is consistent with different observations. 

5.1. Comparison between observed spectral- 
line maps and mass-weighted velocity 
histograms 

The colliding streams we have described as respon- 
sible for cloud formation should manifest themselves 
as multiple peaks in density-weighted velocity his- 
tograms of one velocity component, the equivalent 
in our simulations of the optically thin line profiles 
in observational spectral line maps. Although fre- 
quently Gaussian or other smooth curves are used to 
fit such lines, conferring the idea that the lines are uni- 
modal and correspond to well-defined, isolated enti- 
ties, in actuality the profiles generally exhibit multiple 



peaks (for recent examples, see, e.g., Falgarone et al 
|1994| 



etc. 



which are in fact normally interpreted as 
"clumps" (e.g., WiUiams, Blitz & Stark 1995, here- 
after WBS95). 

For comparison with observational profiles, we show 
in fig. |6| the density-weighted histogram of the 



co mponent of the velocity (ux) for the "parent" cloud 
showed in fig. |l|b. As usual, this histogram has been 
obtained using the local speeds in the cloud frame. 
This histogram may be compared with the ^^CO line 
profile for the Rossete Molecular Cloud (RMC) shown 



in fig. 4 of WBS95^. This cloud is in many aspects 
comparable to the cloud of fig. |^b: while the RMC's 
dimensions (as deduced from fig. 17 in WBS95) are 
- 90 X 70 pc, the cloud in fig.jfj is - 250 x 120 pc. 
Furthermore, the cloud in fig. nb has a mean density 
of ^ 10 cm"'^, while the gas sampled in the spectrum 
of fig. 4 of WBS95 has a mean density ^ 15 cm~^. 

Comparing the histograms for both of these clouds, 
we note several points in common. First, both sets of 
data have FWHMs of roughly 6 km s~^ when only 
the main features are considered. Second, both sets 
exhibit high- velocity bumps, at several km s~^ from 
the centroid. Finally, and probably most importantly, 
the "main" features are seen to contain substructure 
in both sets of plots. However, while such features 
have been traditionally interpreted as "clumps" , in 
our simulations they originate from extended regions 
within the cloud. 

In order to explain this phenomenon, we show in 
fig. |7| a different representation of the parent cloud 
(fig. ^d). The isocontours denote the density field, 
the arrows show the velocity field in the frame of the 
cloud, and the various tones of gray represent zones 
with nearly the same a;-component of the velocity Ux- 
The velocities range from less than —8 km s^^ to more 
than 8 km s~^. Each tone of gray indicates a velocity 
interval of 2 km s~^. Several issues are then observed. 
First, we note that the various peaks in the histogram 
of fig. ^ contain contributions from extended and sep- 
arate regions in space. Second, local density peaks 
often encompass several velocity intervals, in a simi- 
lar fashion to the situation for the entire cloud, sug- 
gesting a kind of self-similarity. Note then that the 
relative height in the histogram depends mainly in the 
mass fraction of the gas with the corresponding value 
of Ux-, since every velocity interval appears to sample a 
wide range of densities. Third, a given feature (peak) 
in the histogram may have contributions from regions 
with very different total velocity vectors, which only 
happen to have a similar a;-component. The above 
results suggest that the identification of "clumps" as 
features in velocity space is risky at best, and possibly 
erroneous, as already pointed out by Adler fc Roberts 



1992 



der 



A final word concernin g th e velocity field is in or- 
As pointed out in § 4T, significant vorticity is 



^Although fig. 4 in WBS95 shows CO spectra, which are opti- 
cally thick, there is little, if any, qualitative difi'erence with the 
optically thin ^'^CO spectra shown in fig. 5 of that paper 
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found within the clouds. Observationally, this may 
be detected as velocity gradients along the cloud's 
projected area on the sky, which indeed have been 
observed (e.g., [Arquilla fc Goldsmith 1986|; iGoodman 
et qd 
& 



1993). Furthermore, Miville-Deschenes, Joncas 



'algarone (1998) have recently reported possible 



direct evidence for vorticity in clouds. 

5.2. Magnetic Field: alignment, advection 
and field reversals 



Concerning the magnetic field topology, the field 
reversals and bendings mentioned in § 4.3 might seem 



contrary to the widely accepted fact that the mag- 
netic field lines projected on the plane of the sky, 
as obtained in polarization analysis (see, e.g., pood- 



tnan et al. 199C), are relatively smooth over the area 



of the clouds. However, on closer inspection of the 
field topology in our simulations, there is no contra- 
diction. As can be seen in fig. ^ the magnetic field 
bendings are correlated with bendings of the density 
features, so that, along elongated density features, 
the field appears rather smooth. Besides, the pro- 
jection along the line of sight can make the field lines 
appear smoother than they a re ( Ballcstcros-Parcdc£ 
md Vazqucz-Scmadcni 1998). On the other hand, 
evidence that the interaction of a shock with dense 
molecular gas is producing a reversal in the line-of- 
sight field component in the Orion/Eridanus region 
has been recently reported (Heiles 1997). 

Note also that it is by now well established that the 
field lines do not in general show a unique direction of 
alignment with the density features ( [Goodman et. al 
1992). In our simulations this is also true. Although 



in particular for the case shown in fig. ^ the field lines 
are mostly parallel to the density features, in other 
cases they are perpendicular (see, e.g., the concave 
region at coordinates x ~ 605 and y ~ 510 in fig. ^). 

5.3. Lifetimes of turbulent clouds 

An important consideration is whether the turbulent- 
fluctuation nature of the clouds proposed here gives 
cloud lifetimes consistent with current estimates. First 
of all, we should point out that clouds are not nec- 
essarily completely disrupted by the turbulence. Fre- 
quently they are simply sheared, split, merged with 
others, or they can just re-expand into the general 
medium if they are not very stronly self-gravitating, 
again indicating that the medium is not in precise 
pressure balance (see the videos mentioned above). 



In any case, a conventional estimate for a typical 
"lifetime" can be given bt r = l/Av, where I is a 
characteristic size of the cloud and Av its internal 
velocity dispersion. As an example, we calculate this 
for the "daughter" cloud of fig. 0c. We find Av = 2.3 
and I ^ 30 pixels — 37.5 pc. Thus, r ~ 



km 

1.6 X 10^ _ 

estimates of GMC lifetimes of order a few xlO'' 



yr. This value is consistent with various 

yr 



(Bash, Green fc Peters 1977; Blitz fc Shu 1980; Larson 



1981; Blitz 1994) 



6. Discussion 

6.1. Applicability to molecular gas and the 
role of 7off 

The analyses from the previous sections have been 
performed on simulations which have reasonably real- 
istic cooling functions and a diffuse heating which de- 
cays with increasing density (mimicking self-shielding) 
( PVSP95 ). As seen from eq. (^), such heating and 
cooling rates imply effective polytropic exponents 
typically smaller than 1. These heating and cooling 
rates are most appropriate for diffuse HI structures, 
which may be very large (up to ~ 1 kpc, although 
small-scale structure must exist there as well). The 
question then arises as to whether the dynamical sce- 
nario we have proposed here is applicable to molec- 
ular clouds (MCs) and their substructure (with sizes 
^ 10 to < 0.1 pc). Although the ultimate test for 
this problem will require new numerical simulations, 
which will be presented elsewhere, we can discuss the 
issue to a certain extent here, based on recent results 



on the role of 7off (Tohhne et al. 1987; VSPP96; Pas 



3ot & Vazquez-Semadeni 1998), which is a measure 
of the thermal characteristics of the flow. 

The main difference between the parameters of the 
simulations presented here and another set more ap- 
propriate for MCs is the effective polytropic expo- 
nent of the flow. (Note that the only scale-dependent 
quantities in our simulations are the heating and cool- 
ing rates. All other quantities are scale-free.) While 
in general in our simulations < 7efi < 0.5 ( c.f. § 
~|), it is possible that in molecular gas 7eff ^ 1 ( |My-| 



ers 1978 ). Actually, the situation at densities ^ 10'^ 
cm~'^ is rather uncertain, with 7off possibly being ~ 1 
between lO'^ and a few xlO'* cm~^ because of the 
importance of coUisional deexcitation and radiative 
trapping in CO lines, and either larger or smaller than 
unity at larger densities depending on the presence of 
embedded stellar sources and other factors (gcalo et 
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al. 1998). 

In any case, we can discuss the expected effect of 
the gas having "f^s ^ 1 on the nature of the den- 
sity fluctuations arising in the flow. It has recently 
been shown by Passot & Vazquez-Semadeni (1998) 
that at larger j^s the density peaks are spatially 
more extended and of smaller amplitude, while the 
density minima are very deep. At small 7cff (typ- 
ically < 7eff < 1), the situation is reversed, the 
peaks being very concentrated and large, while the 
minima are not as deep. This is because the sound 
speed c scales with the density as p(Trff~i)/2^ caus- 
ing the peaks to re-expand faster at larger 7off. At 
7eff = 1 a critical case occurs, the fluctuations of 
Inp bein g symmetrical with respect to the mean. 
VSPP96 have calculated the density jump for shocks 
in a barotropic fluid parameterized by 7cff. Thus, 
if 7cff is larger at MC densities than in our simula- 
tions, we can expect the density peaks to be wider 
and less pronounced. But we do not expect the gen- 
eral scenario of clumps being formed by turbulent 
compressions to change significantly. The colliding 
streams, if supersonic (or super-magnetosonic) , cre- 
ate non-stationary shock-bounded slabs, which sim- 
ply re-expand (if not made self-gravitating by the 
compression-see below) on different timescales de- 
pending on 7off after the compression subsides. If 
the compressions are subsonic no shocks form, but, 
for sufficiently small jcS and/or cloud masses close to 
their Jeans mass, collapse may still be induced by the 
compressions (Tohline et al. 1987). 

One crucial issue, however, is whether the com- 
pressed slabs can become self-gravitating upon the 
compression, and, if so, whether a quasi-hydrostatic 
configuration can be formed. The former problem 



be halted, since Ei increases faster than \Eg\. The 
corresponding values of 7cr for 2- and 1-dimensional 



has been addressed by a number of authors (Euntcr 
197|9|; [Hunter fc Fleck 1982|; |Tohlinc ct al. 1987t Vish- 



niac 1983, 1994; [Elmegreen 1993b| ; |VSPP96D . Two 
main issues are at play here. One is the total gravi- 
tational (Eg) and internal (Ei) energies of the cloud 
(for the purpose of this discussion the "internal" en- 
ergy can be generalized to include all forms of energy 
that provide support for the cloud against gravity). 
The other is the rate at which each energy increases 
upon compression. If \Eg\ > E^, the cloud is gravita- 
tionally unstable, and collapses upon a perturbation. 
This is essen tially an integral version of the Jeans cri- 
terion (e.g., Bonazzola et al. 1987 ). However, for 
three-dimensional collapse, it is well-known that, if 
7off > 7cr = 4/3, then the collapse will ulimately 



collapse are respectively 1 and (VSPP96). 

Conversely, if a cloud is initially stable according 
to Jeans (Ei > l-Bgl), it may be made unstable upon 
compression. Tohline et al. (1987) have calculated 
the required Mach number of a 3-dimensional com- 
pression as a function of 7eff for this to occur, pro- 
vided that 7off < 4/3 (for < 7cff < 1 they found 
that the required Mach number is independent of 
the cloud's mass). Two- and one-dimensional com- 
pressions require 7ofr < 1 and 7cfr < 0, respectively 
(VSPP96). The issue is then the following. In or- 
der for an initially stable cloud to be pushed over the 
internal energy "barrier" (Tohline et al. 1987) by a 



strong enough turbulent compression, it is necessary 
that 7off < 7t, where 7t is the critical 7 for \Eg\ to 
increase faster than E^ upon the turbulent compres- 
sion. From there on, gravitational collapse takes over. 
Since in general turbulent compressions are expected 
to be less than three-dimensional, then generally we 
expect 7t < 7g, where 7g is the critical 7 for gravi- 
tational collapse.^ Then, 7off < 7t < 7g- Thus, we 
conclude that, if collapse was triggered by a turbulent 
compression, it cannot later be halted by the internal 
energy, unless "fcS changes in the process. The lat- 
ter may happen in the late stages of collapse, if the 
dynamical timescales become shorter than the radia- 
tive cooling timescales, causing jcS to approach the 
actual heat capacity ratio of the gas 7, rather than 
the value derived from the condition of thermal equi- 
librium (see, e.g., PVSP9E). 

We should insist here that the effects of all forms 
of support, such as the magnetic pressure and tur- 
bulent support, can in principle be included in jeS- 
On the other hand, this is a very simplistic descrip- 
tion, since in general 7off may be a function of the 
density, and the anisotropic nature of magnetic sup- 
port may not be well described by a simple polytropic 
exponent. In any case, this discussion illustrates our 
main point: quasi-hydrostatic configurations cannot 
be produced out of turbulent fluctuations, unless the 
flow becomes closer to adiabatic than to isothermal 
during the collapse. Fluctuations either rebound or 
colla pse. However, at small 7cff, the rebound time is 
long ( |VSPP96| ; [Passot fc Vazquez-Semadeni 1998[) . 



*Here we are allowing for the possibility of gravitational collapse 
to occur faster in one direction than in the other two, as is the 



case of cosmological "pancakes" 
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We can examine more quantitatively the conditions 
under which adiabatic compression or collapse would 
occur if thermal pressure provides the main support 
mechanism for a clump, by comparing the cooling 
timescale TcooI = nkT/ K (A = cooling rate in units of 
erg cm~'^s~^) to the compression time Tcomp = L/v 
and the collapse time Tcou = (Gp)^^/^. For the 
cooling rate we distinguish two density regimes. For 
10'^ ^ n ^ 10^ cm~^ the cooling is dominated by CO 
line cooling. We adopt a rough approximation for the 



CO coolmg rate (see reterences m bcalo et al. lyys j as 
A ~ iU^^''n3J'^, where ns — n/iU"'cm^'' and we tiave 
assumed that collisional deexcitation yields a linear 
(rather than quadratic) density dependence at these 
densities. We recognize that this is a crude represen- 
tation of the molecular cooling rate, but it suffices for 
our order of magnitude estimates. At larger densities 
the cooling will be dominated by gas-grain collisions 
(if there are no embedded protostellar sources), with 
a rate 2 x lQ-^^r?T^/'^ (T - T^,), where T^, is the 
grain temperature and the coefficient depends some- 
what on the adopted grain parameters. We assume 
T — Tgr K, T for our order of magnitude estimates. 

For compressions, we find that the condition for 
adiabatic evolution (Tcomp < ''cool) yields a critical 
length scale Lpc = L/lpc 



4 X IQ-^M/tH^ 



] 7 X 10-4M/n5 



CO cooling 
gas-grain cooling. 



where M is the Mach number of the compression, 
Tio = T/IOK, and 77,5 — n/10^cm~^. For any rea- 
sonable Mach number adiabatic compression is only 
possible for tiny regions of size < 10^^ pc. 

If a (non-adiabatic) compression results in collapse, 
the condition for adiabatic collapse is t^ow < t^ooI- At 
densities 10'^ ^it- <, 10^, the adopted CO cooling rate 
results in the condition > 2.1 x lO'^rfg- But at 
such large densities the cooling will be controlled by 
gas-grain collisions, not molecular line cooling. Using 
the gas-grain cooling rate, the condition for adiabatic 
collapse is found to be ns < 1.2 x 10~^Tj^^. But 
at such small densities molecular line cooling would 
dominate and be non-adiabatic. Thus for collapse in- 
duced by compression, adiabatic evolution will never 
occur until the collapsing entity is so dense that it 
becomes opaque to its own infrared radiation, which 
only occurs at very large (protostellar) densities. 

These estimates indicate that turbulence-induced 
compressions or collapse will generally be controlled 



by the effective polytropic index jcS set by the ther- 
mal equilibrium condition, since TcooI is smaller than 
the relevant dynamical timescales, except for ex- 
tremely small scales of compression or at protostellar 
densities for collapse. We conclude that quasi-static 
configurations are unlikely over a very wide range of 
physical parameters. 

Is there any observational evidence for a dynamical 
rather than quasistatic state of clumps? A prelimi- 



nary answer appears affirmative. Indeed, Girart et al 



(1997) have recently reported on a system which may 
be interpreted as a collision of streams in L723. Fur- 



thermore, Tafalla et al. (1998) have reported a core 
with inward motions that do not correspond to classi- 
cal gravitational infall models. A detailed comparison 
of these data with our scenario will be presented else- 
where. 

6.2. Effect of dimensionality 

In this paper we have used two-dimensional simu- 
lations, with the known risks involved by them, such 
as a different distance dependence of gravity, the ab- 
sence of a magnetic dynamo, and the possibility of 
reverse cascades (although this is probably less impor- 
tant in the compressible case, where the energy spec- 
trum may be determined by the shocks rather than 
by cascade processes - see, e.g., Vazquez-Scmadcni 
1998| ). This is a necessary price to pay for the larger 



resolution allowed by the two-dimensionality, which 
is essential for the type of hierarchical structure anal- 
yses performed here. 

However, we do not expect the two-dimensionality 
of our simulations to be a matter of concern for the 
validity of our results. The main points raised here 
are those of the mechanism of cloud formation and its 
consequences regarding the mass flux through cloud 
boundaries, the deformation of the magnetic field by 
the flow, and the interpretation of line profile obser- 
vational data. These do not seem to depend on the 
dimensionality of the problem, since they originate 
from the shocks and collisions of streams, which oc- 
cur in both two and three dimensions. In any case, 
intermediate-resolution 3D simulations are currently 
being performed which will allow us to test these ex- 
pectations. 

7. Summary and conclusions 

In this paper we have conducted a close examina- 
tion of the topology of the density, velocity and mag- 
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nctic fields that result in our 2D simulations of the 
ISM at the kpc scale, in particular within the clouds 
that emerge as the turbulent density fluctuations in 
that medium. We found that, in general, the cloud 
boundaries are rather arbitrary, being defined only 
by a density threshold criterion, but not by a physi- 
cal boundary, such as a density or velocity disconti- 
nuity. The same arbitrariness is likely to be present 
in clouds defined by the tracer used to observe them, 
since a cloud may seem to "end" where the tracer 
density falls below a threshold density necessary to 
excite the transition, but this need not correspond to 
any abrupt drop in the mass density. 

We furthermore considered a hierarchy of clouds 
in the simulations, by using a sequence of three pro- 
gressively larger density thresholds for defining them. 
For clouds at all levels of this hierarchy, the veloc- 
ity field in the simulations is continuous across the 
clouds' "boundaries", indicating that the clouds arc 



being formed by colliding gas stre ams ( Hunter et al 



19^; [Elmegreen 1993b| ; |VSPP95| ), with the density 
increasing towards the collision site. The density en- 
hancements need not be confined exclusively to the 
collision site as in Burgers flows, because the veloc- 
ities involved are typically trans- Alfvenic, meaning 
tha|t the information on the coihsion can often prop- 



agate upstream (see also Pouquet, Fassot & Leorat 



19£[T|). Furthermore, since the collisions are generally 
oblique, it is likely that magnetosonic waves prop- 
agating perpendicular to the collision surface may 
overtake the oblique fluid motion, even if the total 
fluid velocity is super- Alfvenic. This scenario is in 
sharp contrast with static models in which the clouds 



arel immersed m an mtercioud medium ot lower den- 
sity but Jiigtier temperature, and contmed by pres- 
sure balance between the two "phas es" (e.g., BM92| ; 
MZ92| ; iMcLaughhn fc Pudritz 199e| ). In our simula- 



tions both "phases" exist, but the transition between 
them is smooth, and the situation, rather than being 
static, is highly dynamic. 

To further characterize the role of the flux across 
the clouds' boundaries, we considered the surface 
(Tkin) and volume (fkin) kinetic terms in the virial 
theorem for the whole ensemble of clouds at one tem- 
poral snapshot in the simulation. It was found that 
both terms are typically of the same order of mag- 
nitude, indicating that the surface terms contribute 
an amount comparable to the total kinetic energy of 
the clouds to their overall virial balance. This was 
interpreted as a signature of the fact that, when the 



clouds are formed by stream collisions, the same flow 
is responsible for both the total kinetic energy con- 
tained within the clouds and for the flux across its 
boundaries. 

We discussed the differences between the thermal 
and turbulent pressures, emphasizing the fact that 
while the former is microscopic and isotropic, the lat- 
ter may involve macroscopic scales, comparable or 
even larger than the clouds' sizes, and is in general 
anisotropic, as indicated by its off-diagonal contri- 
bution to the total pressure tensor Ily . We thus 
concluded that turbulent pressure confinement is a 
contradictory concept, because it involves motions at 
large scales which by definition cause a change in the 
clouds' boundaries. In turn, we argued that thermal 
pressure equilibrium in a non-static medium is irrel- 
evant, because inertial motions may still distort or 
disrupt a cloud, even under thermal pressure balance 
with its surroundings. 

Additionally, we studied the magnetic field topol- 
ogy and strength in the clouds in the simulations, and 
their implication on whether the motions are super- 
or sub- Alfvenic, finding both super- and sub- Alfvenic 
motions within the clouds, supporting the suggestion 



that the field is highly intermittent ( Padoan fc Nord- 
hmd 1998|). We noticed that the field is significantly 



distorted and advected by the fiow, in agreement with 
the result that the flow within the clouds is trans- 
Alfvenic. The field tends to be correlated with the 
density features, although without a unique mode of 
alignment (PVSP95). This causes the field to look 



relatively smooth along density features, in agree- 



ment with most polarization studies (Goodman et al 
1990| ), but at the same time being bent where the 



density filaments are bent as well. 

A crucial question is whether the scenario depicted 
in this paper is consistent with observational data. 
We addressed this problem at the level of density- 
weighted velocity histograms, which can be compared 
to observational spectral line profiles from real clouds 
of similar characteristics to those in the simulations. 
Although there are limitations to this comparison, our 
velocity histograms are encouragingly similar to the 
observational line profiles, both exhibiting compara- 
ble FWHMs (~ 6 km s~^) and high- velocity isolated 
features, and both containing not completely resolved 
substructure, indicating multiple components. We 
emphasized that while this substructure is usually in- 
terpreted as isolated "clumps" , in our simulations it 
originates from different regions of the cloud, which 
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may correspond to different "streams" and in gen- 
eral have one similar component of the velocity but 
very different total velocity vector, as found also by 
Adler & Roberts (1992). We also mentioned obser- 



vational evidence for field reversals in shocks (Heiles 
19!^, and noted that the topology of the magnetic 
field seen in the simulations is consistent with obser- 
vations ( poodman et. al 1992| ). 

Finally, we discussed at considerable length the 
possibility that the scenario of clouds as turbulent 
density fluctuations may apply to molecular clouds 
and their substructure (clumps and cores). Essen- 
tially, only the adopted cooling and heating functions 
in our simulations distinguish them from a case more 
appropriate to molecular clouds, which are believed 
to behave in a nearly isothermal manner. Based on 
recent results on the effect of the effective thermal 



be havior ([rohhne et al 1987| ; |VSPP96| ; [Passot fc 



Va: :quez-Scmadcni 1998| ), parameterized by the ef- 
fective polytropic exponent 7off, we conjectured that 
only the shape (height and width) of the density fluc- 
tuations is expected to change, but not the general 
character of the inflowing streams characteristic of 
their turbulent nature. We also suggested that quasi- 
hydrostatic configurations cannot be produced from 
turbulent compressions, and that density fluctuations 
must either collapse or rebound, unless the effective 
thermal behavior becomes closer to adiabatic during 
the collapse (7off > 4/3). However, for small 7off, the 
rebound times can be quite long. We demonstrated, 
using appropriate cooling rates, that adiabatic evo- 
lution, and hence quasistatic configurations, are un- 
likely to occur either during compressions or during 
compression-induced collapse. 

In summary, we conclude this paper by stating that 
the results presented here are not in contradiction, 
but rather quite consistent, with observational data. 
They seem to be, however, discrepant with the stan- 
dard interpretation of those data in terms of sharply- 
bounded clouds in pressure equilibrium with their sur- 
roundings, and with the identification of velocity fea- 
tures in spectral-line data with isolated "clumps" . It 
is interesting that our results give theoretical support 



to the suggestion made long ago by Chandrasekhai 



& . vlunch (1952^ that the interstellar medium may 
not best be visualized as a system of discrete clouds. 
As they pointed out, the "cloud model" is difficult 
to displace because, as they say, "...a tendency to ar- 
gue in circles can be noted in the literature, in that 
confirmation for the picture of interstellar matter as 



occurring in the form of discrete clouds is sought in 
the data analyzed" . 
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Fig. 1. — Three hierarchical levels of clouds in the 
simulation, showing the logarithmic density scale with 
shades of gray, and the velocity field with arrows. 
Pixel numbers are indicated by the grids, a) The 
whole simulation field field, 800 pixels on a side, cor- 
responding to 1 kpc. The gray scale shows structures 
with p > pt = 1-5 cm""^ b) Magnification of the com- 
plex that appears when setting = 4 in the upper 
right quadrant of fig. la. c) Magnification of the cloud 
that appears when setting = 8 in the lower left re- 
gion of fig. lb. Black arrows indicate super- Alfvenic 
velocities in the frame moving with the cloud, and 
white arrows indicate sub- Alfvenic velocities. The 
contours give the magnetic field strength B. Thicker 
countours indicate larger B-values. The maximum 
density value is 55 cm~^. 
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Fig. 2. — Cuts passing through the cloud shown in 
fig. |l|c, (as indicated by the Unas in fig. ||), showing 
the log of the density (solid line) , and the x- (dotted 
line) and y- (dashed line) components of the velocity, 
a) Cut along x at y — 514. b) Cut along y a.t x = 
602. c) Cut along y at x — 614. In a), two peaks 
are indicated. Peak 1 is seen to correspond to an 
x-shock (abrupt negative gradient) in a) and to a y- 
shock in b). Peak 2 is almost imperceptible in u^, but 
is seen to correspond to a y-shock in c). Also, note 
the "plateau" in Ux extending 50 pixels to the right of 
peak 1, with average negative values, although with 
further substructure. Note that the scale is larger 
than in fig. in order to give a broader view of the 
velocity components. 
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Fig. 3. — The same cloud as in fig. 0c, but with super- 
imposed velocity divergence contours. Density max- 
ima are seen to correspond to maximum negative val- 
ues of the divergence. Also indicated are the lines 
along which the cuts of fig. || are taken. 
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Fig. 4. — Same cloud as in fig. |l|c, but with superim- 
posed vorticity contours. There is a significant corre- 
lation between the density (and therefore the diver- 
gence) and the vorticity fields. We interpret this as 
a result of the generally oblique nature of the stream 
collisions. Solid (dashed) contours indicate positive 
(negative) vorticity. Also shown are the magnetic 
field vectors (arrows). Note the correlation between 
density and magnetic field bendings and the field re- 
versal along the central ridge of the cloud. Note also 
the general alignment of the field with the density, 
although at many places the field crosses the cloud's 
"boundary" perpendicularly (e.g., at the concave re- 
gion below the left filament). 
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Fig. 5. — Kinetic energy £'kin vs. the kinetic surface 
term Tkin for all clouds with areas larger than 300 
pixels in the whole field. The complex of fig. 0b is 
shown by a circle, and the cloud in fig. |l|c is shown 
by the star. Note the similarity of the two terms for 
most of the clouds, indicating that both contribute in 
similar amounts to the clouds' virial balance. 
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Fig. 7. — Same cloud as in fig. |l|b, but indicating 
the magnitude of the local a;-component of the ve- 
locity (gray-scale), in order to identify the origin of 
the x-velocity features appearing in the histogram of 
fig. ^. The velocity ranges are indicated in the side- 
bar. The arrows and contours indicate the velocity 
and the density fields, respectively. Three points are 
noticed. 1. The contribution to any given velocity 
interval originates from extended regions throughout 
the cloud. 2. Regions contributing to a given velocity 
interval contain zones of significantly different densi- 
ties. 3. Regions with the same x-component of the 
velocity can have completely different total velocity 
vectors. 



Fig. 6. — Density-weighted velocity histograms for 
the parent cloud (fig. |l|b). Note the FWHM -5-7 
km s~^, and the multi-component structure. Also, 
note the low-decaying wings. 
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